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W)' A differential geometric approach to singular perturbation theory is pre- 

' sented. It is shown that singular perturbation problems such as multiple-scale 

and boundary layer problems can be treated more easily on a differential ge- 
p) • ometric basis. A general method is proposed based on differential forms and 



Lie-derivatives. Examples from multiple scale theory, boundary layer theory 



^ ' and WKB-theory are given and it is demonstrated that without the a priori 



knowledge of the scaling behaviour of the problem the correct asymptotic ex- 
pansion can be derived with the aid of differential forms. The method is well 

o 

' suited for a mechanical implementation in computer algebra programs. 
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^ . I. INTRODUCTION 

Ph. 

> : 

^ , The analysis of boundary layer problems and multiple scale phenomena which are gen- 

H : 

5^ , eralized under the notion of singular perturbation problems has played a significant role 
in applied mathematics and theoretical physics [1,2]. Regular perturbation theory is often 
not applicable to various problems due to resonance effects or the cancellation of degrees of 
freedom. In order to obtain a uniformly valid asymptotic expansion of a solution for these 
singular problems a whole bunch of methods has been developed such as boundary layer 
expansions, multiple scale methods, asymptotic matching, stretched coordinates, averaging 
and WKB-expansions. Although these methods work well for the respective general 
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theory that unifies all these methods is still lacking [4] . A general applicable method would 
be highly desirable in order to obtain a mechanical scheme that can be easily implemented 
in algorithmic programming languages. 

A physical system often involves multiple temporal or spatial scales on which character- 
istics of the system change. In some cases the long time behaviour of the system can depend 
on slowly changing time scales which have to be identified in order to apply multiple scale 
theory. The choice of the slowly or fast changing scales is a nontrivial task, which presumes 
a good understanding of the physical behaviour of the system and can sometimes only be 
justified by the final result. A naive expansion in a power series of the small parameter is 
often prevented by the appearance of resonant terms in higher orders. These terms have to 
be compensated by the introduction of counterterms. 

Boundary layers are also a common feature of singular perturbed systems. In these 
cases higher order derivatives disappear in the unperturbed equations which leads to the 
cancellation of degrees of freedom of the system and finally in small regions where the system 
changes rapidly. The identification of the fast changing scales is crucial for the solution of 
the perturbation problem and is a subtle task, sometimes involving delicate arguments. 

In this Letter a geometrical approach is presented which allows a mechanical algorithm for 
the determination of the rapid or slow scales on which the system changes. For that purpose 
the system of ordinary differential equations is interpreted as a system of differential forms. 
It is shown that a small perturbation of the system of differential forms can be viewed as a 
deformation of the solution manifold in configuration space. By searching for a Lie-derivative 
that transforms the unperturbed differential forms into the perturbed forms a mapping in 
configuration space is defined that can be used to map the solution of the unperturbed system 
into the an approximate solution of the perturbed system in order to obtain uniformly valid 
asymptotic expansions of the solution and thus determining the secular scales. The method 
is exemplified for common problems like multiple scale problems, boundary layer problems 
and WKB-expansions. 
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II. GENERAL THEORY 



Assume a system of ordinary differential equations in tfie n-dimensional configuration 
space i?" of the following form: 

where the dot denotes the derivative by the independent parameter t. The system is 
solved by a n-dimensional curve C : Xi{t) parametrized by t. It is easily to shown that all 
ordinary differential equations can be written in this form [3,5]. The system of equations 
can be written as a system of equations for differential forms: 

OUij = Fi {Xi, Xn) dXj - Fj {Xi, ...,Xn) dxi ^ (2) 

The differential forms Uij establish a basis for a subspace S of the space of all differential 
forms in i?". Every solution curve C then has the following property: 

jf{x)uj = Q (3) 
c 

for all elements uj of the subspace S and for all scalar functions / on i?". This property 
can be used as a very general definition of a solution manifold of a system of differential 
forms and can be easily generalized to higher order differential forms and closed ideals of 
differential forms. 

Now we consider a perturbed subspace S' of a subspace 5*0. Let be a basis of an 
unperturbed subspace 5*0 and o;*^*^ = uJq^ + eujf' be a basis of the perturbed subspace S' 
where £ is a small parameter. The condition for a solution curve of the perturbed system S' 
can be written as: 

j f{x){u;o + eu;,)^0 (4) 
c 

for all elements lo — loq + eui of the subspace S' and for all scalar functions / on i?". 
We assume that there exists a mapping $ of i?" into i?" that maps the solution curve Co of 



3 



the unperturbed system into the solution curve C — ^ (Co) of the perturbed system. If the 
mapping fulfills the condition: 

(cuf ) + ecoP) = uip + 6 J2 ^^'M'^ (5) 

j 

with appropriate functions A^*'''^ (x) then one obtains: 

J f = / (/a;«) = E / ('^'* (/) + (/)) 4'^ = (6) 

*(Co) Co ^' Co 

By expanding the mapping in a power series with the small parameter e one obtains 
— l + eL + O (£^) where L is an appropriate Lie-derivative, one obtains in first order: 

u;P + l4^^^X('M'^ (7) 
j 

This is a linear equation for the Lie-derivative L which we can try to solve. It defines Lie- 
derivatives that map solutions of the unperturbed system 5*0 into solutions of the perturbed 
system S'. This method can be applied also to singular perturbation problems where a 
straightforward regular expansion is not available. 

III. BOUNDARY LAYER ANALYSIS 

The first example of a singular problem that can be solved with the new method is given 
by the following boundary layer problem [1] : 

ey" + y' + y^O (8) 

with £ ^ 1 and the boundary conditions y (0) = and y' (0) = 1 . The prime denotes 
differentiation by x. This is a common boundary layer problem introducing a layer of 
thickness e. The layer thickness is usually obtained from a dominant balance argument. 
The system is then rescaled to obtain a inner region equation which is afterwards matched 
to the outer region solution. The choice of the thickness of the boundary layer is justified 
in the end of the calculation by the existence of a uniform solution. 
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In the proposed method such a rescahng is not necessary.. The original equation is 
written as a system of differential forms: 

dy — z dx = (9) 



dy + y dx = —e dz 



(10) 



The unperturbed differential forms are given hy ujq = dy — z dx, ujq = dy + y dx and 
uj[^'^ = 0, uj^'^ = —dz are the perturbation, respectively. The unperturbed system can be 
rewritten as dx = and dy = which has the solution x = xo and y = yo. This zero-order 
solution is clearly a singular solution. In deviates strongly from the zero order solution 
usually obtained from the equation y' + y = 0. This difference is crucial in the present 
approach. 

After a tedious but straightforward calculation one obtains a Lie-derivative that maps 
the unperturbed differential forms into the perturbation: 



L = \n {y + z) + {z - y\n {y + z)) d. 



y 



(11) 



One easily checks the transformation properties: 



(1 + eL) {dy + y dx) = dy + y dx + e dz + efiUii 



(0) 



(12) 



(1 + eL) {dy — z dx) = dy — z dx + egiOJl 



(0) 



(13) 



where /j and Qi are functions of x, y and z as defined in equation (^. From the zero-order 



solution X = Xq one obtains: 



(1 + eL) (x — Xq) = X — Xq + e\n.{ii + z) 



(14) 



After the substitution z = y' the solution can be written as: 



X = Xq - e\n{y + y') 



(15) 



5 



and thus by substituting and A = exp (xq/e) one obtains the following first-order differ- 
ential equation: 

y' + y = Ae"^/^ (16) 

The effect of the mapping from the unperturbed system to the perturbed system is now 
clearly demonstrated. In the limit e —>■ the term with the second derivative in equation 
(P) vanishes and the order of the equation changes which implicates a boundary layer which 
becomes infinitely thin for small e. The transformed system (|1^) still possesses the additional 
degree of freedom for small e and one observes how the singularity forms when e becomes 
small. Equation ([T6|) can be easily solved in order to obtain 

y = Ae"^ J e^(^-^/")rfx + Be-"" (17) 

and finally 

y = Ae-^'l^ + 5e-^ (18) 

which is the correct general uniform first-order approximation. It should be emphasized 
that the present approach gives the correct scaling behaviour in the boundary layer without 
the need for a rescaling of the independent variable which is necessary in the usual asymptotic 
matching procedure or in the renormalization group approach. Thus the presented method 
is superior to the conventional methods in that it chooses selfconsistently the scaling of the 
independent and dependent variables in order to obtain a regular expansion. 



IV. MULTIPLE SCALES 

Another problem very common in mathematical physics is the appearance of multiple 
scales on which the system changes. An example for this behaviour is given by the following 
linear oscillator with a nonlinear damping term [1]: 

y"^y^ey"' = {) (19) 
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The prime shall denote differentiation by the time variable t . The system possesses a 
short time scale of the oscillation and a long time scale of the damping. The amplitude of 
the oscillation slowly decays with time and a frequency shift might be observed. 

The unperturbed system is just an ordinary linear oscillator with frequency unity. The 
zero-order solution is given hj y = R cos (t + 6) with arbitrary constants R and 6. In order 
to simplify the calculations we introduce new dependent variables R and 6 defined hj y = 
R cos (t + 9) and y' = —R sin {t + 9) . A basis for the perturbed subspace S' is given by: 

dR - eR^ (--- cos (2 (t + 9)) + - cos (4 {t + 9))] dt (20) 
\8 2 8 / 

d9 + eR^ Q sin (2 (t + 0)) - ^ sin (4 {t + 9))^ dt (21) 

The zero-order forms are just dR and d9 while the perturbations are given by the terms 
with the front factor e of equations (pO|) and (|2l|). The system can further be simplified by 
introducing the variable u = 1/R^ giving rise to the basis of the perturbed system: 

du-2e{\-\ cos (2 {t + 9)) + \ cos (4 (t + 9))] dt (22) 
\8 2 8 / 

ud9 + sin (2 (t + 9))-^ sin (4 (t + 9))^ dt (23) 

A Lie-derivative that transforms the zero-order differential forms into the perturbing 
differential forms can be obtained as: 

- i sin (2 (t + 0)) + ^ sin (4 {t + 9))^ + (24) 
i Q cos (2 (t + 0))-^ cos (4 (t + 9))^ de (25) 

The transformed zero-order solution is then given by: 

u = UQ + e{h-]^ sin (2 (t + 0)) + ^ sin (4 (t + 9))j (26) 

= ^^0 - ^- cos (2 (t + ^^)) - 4 cos (4 (t + 9))) (27) 
uo \8 61 J 



For large t the sin and cos terms can be neglected and one obtains as the final approxi- 
mation: 

, , (28) 

which is the correct uniform approximation for large times as obtained from the multiple 
scale method. The time scale et enters naturally by the differential forms method and there 
is no need for an a priori rescaling of the time variable or the introduction of multiple time 
scales. 



V. WKB-EXPANSION 

The last example we give is also a very standard problem in mathematical physics and 
one of the building blocks of wave mechanics, the WKB-approximation of linear differen- 
tial equations. The standard WKB-problem is given by the following ordinary differential 
equation [1]: 

eY + ^''{x)y = Q (29) 

The basis of differential forms of the unperturbed subspace 5*0 associated to this problem 
and the perturbing forms are given by — dy — z dx, Uq^ — (x) y dx and — 0, 
oji = dz. The basis can be rewritten in order to exhibit the exact nature of the zero-order 
forms as: 

u;W=dx + e'-^ (30) 
Q2 (^x) y 

"™ = ''(''Va4)<-') (31) 

By introducing the new variables u — y^ and v — z^ a Lie-derivative that transforms the 
zero-order basis into the perturbation can be found which reads as: 

1 ^ ^d, + v^^ (32) 



n^{x) VV u 



By using the zero-order solutions x — and y — one obtains for the first-order per- 
turbed equations: 

y'^±'-n{x)y (33) 
and thus recovers the WKB-approximation for the initial problem: 

^ = e±i/"(^)'^^ (34) 

This solution shows again that the differential geometric method simplifies the task of 
obtaining a uniform solution for a singular problem. 



VI. CONCLUSIONS 

In this Letter we have shown that the differential geometric method is a superior tool in 
order to obtain uniformly valid asymptotic expansions of singular perturbation problems. A 
basic problem of singular perturbation theory is the appearance of different scales on which 
the dependent variables change. The identification of these scales is a nontrivial process 
which requires a deep understanding of the problem under consideration. The correct choice 
of the scaling of the independent variables strongly influences the region of validity of the 
asymptotic solution. The appearance of multiple scales can further complicate the task 
of finding a uniformly valid approximation of the solution. The method presented in this 
Letter facilitates this task and gives the correct expansions in a very natural way. We 
have given three classical examples from singular perturbation theory and have shown how 
easily uniformly valid approximations can be obtained with the new method. The method 
is based on very general geometrical considerations and thus can be easily extended to more 
complicated problems involving partial differential equations (work in progress). 
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